Chromatin dynamics: Nucleosomes go mobile through twist defects 
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We study the spontaneous "sliding" of histone spools (nucleosomes) along DNA as a result of 
thermally activated single base pair twist defects. To this end we map the system onto a suitably 
extended Frenkel-Kontorova model. Combining results from several recent experiments we are able 
to estimate the nucleosome mobility without adjustable parameters. Our model shows also how the 
local mobility is intimately linked to the underlying base pair sequence. 

PACS numbers: 87.15.He, 36.20.Ey 



The genetic information of all higher organisms is orga- 
nized in huge beads-on-a-chain arrays consisting of cen- 
timeters to meters of DNA wrapped around globular ag- 
gregates of so-called histone proteins. The basic unit 
of chromatin, the nucleosome, is a tiny 10 x 5 x 6 nm 
sized spool composed of 147 base pairs (bps) DNA 
tightly wrapped around an octamer made from 8 histone 
monomers. Each nucleosome is connected via a stretch 
of "linker" DNA to the next such protein spool. The 
wrapped DNA, being coiled in ~ l| turns of a left handed 
helix with radius ~ 4.2 nm, is strongly distorted from 
its preferred straight ground state due to strong interac- 
tions with the histone octamer, namely short range elec- 
trostatics (between the negatively charged DNA sugar- 
phosphate backbone and the positively charged octamer 
surface) and through extensive hydrogen-bonding - both 
localized at 14 discrete interactionpatches helically ar- 
ranged along the octamer surface [jj]. 

Higher order structures, from the 30nm-chromatin 
fiber up to the highest level of DNA condensation, the 
fully folded chromosome, are designed to achieve a huge 
DNA volume fraction. They all rely on the signifi- 
cant stability of the nucleosome complex. On the other 
hand, fundamental life processes like transcription (mak- 
ing RNA offprints from the underlying DNA) and DNA 
replication seem to be in conflict with the picture of a 
stable nucleosome, as they are all performed by protein 
machines that track the DNA helix. The latter inevitably 
implies that every DNA bound obstacle (protein) has to 
be penetrated or even completely removed from its DNA 
target. In fact, the numbers are quite dramatic: A typi- 
cal gene extends over hundreds of nucleosomes, each con- 
tributing 30 — 40 ksT net adsorption energy Also 
other mechanisms like the activation of genes rely on reg- 
ulatory protein binding to specific DNA sequences that 
are often covered by nucleosomes making them inacces- 
sible. 

A key to the understanding of these seemingly con- 
tradictory features might be the physical phenomenon of 
thermally driven nucleosome "sliding" along DNA (also 
called nucleosome repositioning) which has repeatedly 
observed in well-defined in vitro experiments 0, El E| , re- 
viewed in Ref . Q ■ Spontaneous repositioning is strongly 



temperature dependent; at room temperature nucleo- 
somes move a few tens of bps within an hour || . Despite 
clear evidence for repositioning the underlying mecha- 
nism has been the matter of longstanding controversy, 
especially due to the lack of any quantitative theoreti- 
cal treatment of nucleosome statics and dynamics that 
has to rely on the detailed knowledge of the molecular 
structure and its underlying parameters. 

Only very recently - since the documentation of the 
high resolution X-ray structure_ll| and the presentation 
of other new experiments 0, 0, @ - this has become pos- 
sible. First theoretical models of nucleosome reposition- 
ing 0- 113 assume that it is based on the formation of 
DNA loop defects that form on either end of the nucle- 
osomal DNA followed by their thermal diffusion around 
the octamer, similar to the de Gennes-Edwards reptation 
mechanism. This model seems to be successful in explain- 
ing the apparent 10 bps quantization of the nucleosome 
"jump" length || and it also reproduces the observed 
diffusion constants. Recent experiments 0, however, in- 
dicate a more local 1 bp-step mechanism that cannot be 
understood within this model. This lead us here to con- 
sider an alternative mechanism: twist diffusion. The car- 
rier of motion in this case is a twist defect that contains 
one missing or one extra bp. 

The possibility of twist defects was demonstrated as 
soon as the high resolution crystal structure of the 
core particle (the octamer plus wrapped DNA) was re- 
solved 0. In that study the core particles were recon- 
stituted from palindromic 146 bp DNA and histones as- 
suming that this would result in a complex with perfect 
two-fold symmetry. However, it turned out that one bp 
is localized at the dyad axis, the rest being divided into a 
73 bp half and a 72 bp half. The missing bp of the shorter 
half is, however, not localized at its terminus but instead 
at a 10 bp stretch close to the dyad axis (cf. Fig. 4d in 
Ref. 0). This is due to the attraction between the DNA 
termini of adjacent particles in the crystal that come close 
to mimic a bp step at the cost of forming a twist defect 
far inside the wrapped chain portion. This allows us to 
estimate the energy for a single defect to be smaller than 
the stacking energy that is ~ 10 - 20fc B T Q. 

In order to model the twist diffusion mechanism we 
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map the nucleosomal DNA on a Frenkel-Kontorova (FK) 
chain of particles connected by harmonic springs in a 
spatially periodic potential (cf. Fig. 1). The original FK 
model was introduced more than sixty year s ag o to de- 
scribe the motion of dislocations in crystals [llj . In the 
meantime variants of this model were applied to many 
different problems including charge density waves |12j . 
sliding friction 0, Q] , ionic conductors 0, 0] , chains 
of coupled Josephson junctions and adsorbed atomic 
monolayers [ill llflj . Here, in the context of DNA ad- 
sorbed on the octamer, the beads represent the base 
pairs. The springs in between have an equilibrium dis- 
tance b — 0.34 nm and a constant that reflects the cou- 
pled DNA twist-stretch elasticity. Specifically 

Elastic ({Xn}) =Y,C [ ^ k+1 ~ ^ l) 2 (1) 

Here the conformation of the wrapped DNA is given by 
the set {x n } where x n is the position of the nth bp mea- 
sured along the helical backbone; C ~ 70 — lOOfcsT is 
the combined twist and stretch spring constant including 
the (here unfavorable) twist-stretch coupling [2fjj and the 
summation goes over all bp associated with the wrapped 
DNA. In addition there is the external potential of the 
14 contact points to the octamer with neighboring points 
being 10 bp apart [l| that we model as follows 



Eads ({x n }) = - 




x9 (a - \x k - Wbl\) (2) 

with 8 being the Heaviside step function. The two pa- 
rameters of the external potential, its depth Uo and its 
width a, can be estimated as follows. Uq represents 
the pure adsorption energy per point contact which fol- 
lows from competitive protein binding [2] to be of or- 
der 6/csT. The other parameter, a, can be estimated 
from the fluctuations of the DNA in the crystal mea- 
sured by the B-factor (cf. Fig. lb in Q]) at different mi- 
cleosome positions. The ratio of DNA helix fluctuations 
Rfiuct = (x 2 middle ) I (x 2 bond ) « 3 at positions between 
the binding sites and at the bound sites is a measure of 
DNA localization. Using a quadratic expansion of Eq. El 
one finds from a straightforward normal mode analysis 
that a = (5U /[(Rfiuct ~ l)C]) 1/2 6 - b/2, i.e., the ad- 
sorption regions lead to a strong localization of the DNA. 
Knowing all involved parameters the total energy of the 
DNA chain confined in the nucleosome can be written 
down 

Etot — E e i as tic + E ads + E sd (3) 

The last term E sd is the sequence dependent part of the 
total energy which we will neglect first. In the follow- 
ing we study the mechanism for thermal motion of DNA 




FIG. 1: The twist-diffusion mechanism for nucleosome repo- 
sitioning, a) A concerted translational and rotational motion 
of DNA leads to injection of twist-defects (kinks) which mi- 
grate between the octamer adsorption sites (black triangles) 
leading to a "creep" motion of DNA. b) The corresponding 
Frenkel-Kontorova model for twist diffusion and its charac- 
teristic parameters (cf. text for details). 



governed by E to t- Generally two scenarios are possible: 
(?) The generation of kink-antikink pairs inside the nu- 
cleosome and (tt) a kink (or antikink) injection at either 
nucleosome end. Since the first mechanism is energeti- 
cally roughly twice as costly than the second one, we will 
focus here on the (anti)kink injection mechanism only. 

How and how fast does the kink step around the nu- 
cleosome? Due to the strong DNA localization at the 
binding sites (a/b < 1) for a realistic range of parame- 
ters Uo and C the kink is localized either between two 
adsorption positions, i.e., smeared out over 10 bp (de- 
noted by the K w state), or between three of them, i.e., 
smeared out over 20 bp (the K20 state). It is obvious that 
the motion of a (anti)kink will consist of an alternation 
between K\o and K20 states similarly to an earthworm 
creep motion. To model this process we introduce the 
effective kink coordinate xk describing the coordinate of 
the DNA bp being pinncd/depinned during a single kink 
step, so that xk ~ and Xk ~ b/2 correspond to K w 
and K20, respectively, whereas xk ~ b means that the 
kink moved by one bp step. The Peierls-Nabarro poten- 
tial experienced by the kink is then given by Upn (%k) = 
C eff (x K /b - 1/2) 2 - U (x 2 K /a 2 - if for < x K < b/2 
and Upn (xk) — Upn {b — xk) for b/2 < xk < b. Here 
Ceff = io±i C w hh the "— " sign referring to a kink (1 
bp missing) and the "+"sign to an antikink (1 additional 
bp). Depending on the ratio of parameters Uo and C, 
the state if 20 corresponds to a local minimum or max- 
imum of Upn whereas K 10 is always stable for the rel- 
evant parameter range. The rate for the kink step pro- 
cess is then given by the expression fstep = kp,Tjo/b 2 ( e ff 

With J' 1 = (jj e -UPN(sb)/k B T ds ^ e +U PN (sb)/k B T d ^ 

and ( e ff = j^Hspin, the effective kink friction constant. 
Here fj, sp i n = 1.3 x 10~ 20 Ns is roughly the rotational fric- 
tion for a single basestep |2l|. To determine the rate at 
which twist defects are formed at the entry/exit points of 
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FIG. 2: The "charged" Frenkel-Kontorova extension of the 
model in Fig. lb) including effects from anisotropic sequences. 

a) In addition to the contact points the DNA sequence cou- 
ples linearly to an octamer- fixed "bending field" through 
the anisotropic bending parameters qi ("bending charge"). 

b) Two sequences with extremely different mobilities. SI: 
highly anisotropic, 10 bp phased ("TG"-like) sequence with 
D sc i ~ 10~ 4 — 10~ 5 bp 2 /s. S2: random sequence correspond- 
ing to > 95% of the genome with D a d ~ 10 2 bp 2 /s. 



the DNA one can now use an argument similar to the one 
presented in Ref. : The ratio of the life time tu / e of a 
kink to the time interval t\ n j between two kink injection 
events at the end of the wrapped DNA portion equals 
the probability to find a defect on the nucleosome, i.e. 



tlife/t 



in j 



N sited 



-Uu 



Here N site = 13 denotes 



the number of possible positions of the defect between 
the 14 binding sites and Uxink — C/10 is the energetic 
cost for a single kink (cf. above). 

How is the average life time tuf e of a defect related 
to t step , the typical time needed for one step? This can 
be determined from the mean first passage times r/ e / t 
and T r ight for a defect that forms, say, at the left end to 
leave the nucleosome at the same or at the other end, 
respectively. From Ref. one finds 7} e /t = (25/6)£ S i ep 
and T r ight = 28t step - Furthermore, the probability to 
leave at the left end is V i e ft = 12/13 and at the right 
end Pright = 1/13 |22j which gives the life time as 
the weighted average tuf e = 6t ste p- Only a fraction 
Pright of the defects reaches the other end and will lead 
to a repositioning step, i.e., the time of a lbp diffu- 
sion step of the nucleosome along the DNA is given by 
T = ti n j / Pright- Putting all this together we arrive at 
T ~ 6& 2 Ce#J,7 1 / fc B' r exp(C/10fc B T). For realistic pa- 
rameter values C — lOOfceT, Uq = 6fceT and Rfi uc t — 3 
we find T ~ 10~ 3 s implying a nucleosome diffusion con- 
stant D = 580 bp 2 /s = 6.6 x 1CT 17 m 2 /s. 

Hence we find repositioning rates that are orders 
of magnitude faster than the ones observed in experi- 



ments 0. Even worse, the experimental observation of 
an apparent 10 bp jump length [j| seems to be inconsis- 
tent with our predictions. We show now how these facts 
can be explained by the existence of additional barriers 
with a 10 bp periodicity. To do so we have to extend 
our simple model to deal with the quenched disorder 
stored in the DNA bp sequence. The sequence depen- 
dent anisotropic bendability, i.e., the propensity of DNA 
to bend in different directions with different elastic con- 
stants turns out to be essential. It has been known for 
long HI that (A/T) rich dinucleotide steps (dns) 
prefer to face the octamer in the minor groove (i.e., at 
the octamer contact points) whereas (G/C) rich dns pre- 
fer to face the octamer in the major groove (i.e., between 
contact points) . This reflects different propensities of the 
dinucleotides to widen or compress towards the DNA mi- 
nor groove. To incorporate these anisotropic effects into 
our model we first note that the bending state of the DNA 
molecule is fully constrained by its helical path on the 
octamer surface. Moving a DNA sequence via twist dif- 
fusion by a few bp (< 10 bp) along that path changes the 
relative rotational setting of the bent DNA with respect 
to its preferred bending direction causing an energetic 
penalty, whereas a motion by 10 bp restores the initial 
rotational setting. We address this by introducing a lObp 
periodic "bending field" Fb en d (x) = — cos [2nx/ (10b)] 
attached to the octamer surface. We assume the DNA 
sequence to couple linearly to that field through "bend- 
ing charges" qk attached to each of the dns. This gives 
us finally the third term in Eq. 



E s d — ^ qkFbend (Xk) + ™fe 



(4) 



In addition to the anisotropic term we also introduced 
here the isotropic bending parameters rrik to include 
isotropic flexibility effects (which become important 
when the qk 's vanish or average out). The summation 
involved is again over all base pairs |25j incorporated in 
the nucleosome. qk and rrik both have units of energy and 
can be extracted from competitive protein binding exper- 
iments for each of the 10 dns (AA, AT, GC...). To 
obtain a rough estimate we distribute the dns into three 
classes: 1) (G/C) containing dns, 2) (A/T) containing 
dns and 3) mixed dns (like AG, CT etc.) and treat the 
dns in each class as identical. Using the available ex- 
perimental data [3 Ei| we then arrive at qc/c ~ 95, 
q A/T w -85 , q mlxe d ~ and m G/c w 20, m A/T « -3, 
m-mixed ~ 7, where all energies are in cal/mol per dns. 

It turns out that the nucleosome mobility depends 
strongly on the underlying bp sequence. When shifting 
the position of all beads by I bp steps, Xk — ► Xk + lb, we 
find E s d (I) — (A/2) cos (27rZ/10 — <j>) to vary as a cosine 
function of I with phase <j) and amplitude A determined 
by the DNA sequence, which is assumed to be appro- 
priately periodic here. Arranging G/C and A/T tracts 
properly and taking the sequence dependent q and m 
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values given above we can easily reach amplitudes A (i.e. 
barriers to repositioning) that exceed 10 — 12 kcal/mol! 
A very effective sequence arrangement called the "TG"- 
sequence which leads to a strong nuclcosome stability and 
localization was experimentally constructed in Ref. 
by putting G/C tracts around positions k = 0,10,20... 
and A/T tracts around k = 5,15,25.... In our picture 
this means to put the "bending charges" q along the 
DNA such that they couple favorably to the bending field 
Fbend for a distinct rotational setting whereas a 5 bp shift 
is extremely costly (cf. Fig. 2). The 5S-RNA sequence 
which was used in most nucleosome mobility experiments 
shows also the effect of an optimal rotational setting. It 
is less pronounced than in the "TG" case, yet it is still de- 
tectable. More involved theoretical computations relying 
on molecular sequence dependent deformability parame- 
ters [27| reveal barriers A w 5 — 6 kcal/mol for this par- 
ticular sequence. The sequence dependent barrier height 
A exponentially suppresses the bare (sequence indepen- 
dent) diffusion constant D obtained above leading to the 
sequence dependent diffusion constant D s d'- 

D sd = DI ( f (A/2k B T) « I£± e ~(A + c/io )/kB T (5) 

with Iq being the modified Bessel function. 

Equation predicts that mobility experiments with 
highly anisotropic sequences like "TG" (instead of the 
standard "5S-RNA") would find hardly any appreciable 
repositioning on the one hour timescale if it would be 
solely mediated via twist defects (D sc i = 10~ 6 — 10~ 7 x 
D = 10- 4 - 10- 5 bp 2 /s). The typical path for a nu- 
cleosome to escape from such a rotational trap goes very 
likely via the previously considered loop formation mech- 
anism that allows "tunneling" over sequence bar- 
riers, thus dominating over twist-diffusion for extremely 
anisotropic sequences. An experimental test for this pre- 
diction would be to increase the free DNA segment length 
which in this regime should strongly enhance the loop 
mediated mobility [l(i| whereas it would leave the twist 
diffusion unaffected. Going to the other extreme, in the 
most relevant case of random isotropically bendable se- 
quences which make up more than 95% of the eucaryotic 
genome one should observe that the twist diffusion mech- 
anism is strongly enhanced by 2-3 orders of magnitude 
as compared to the in vitro measurements on " 5S-RNA" . 

In conclusion the following picture is implied: On phys- 
iological timescales the majority of genomic nucleosomes 
seems to be intrinsically highly mobile. However, only a 
small fraction (< 5%) of all nucleosomes has strongly re- 
duced mobility due to anisotropic DNA sequences which 
they populate. We speculate that only the latter re- 
quire the action of active (ATP consuming) remodelling 
mechanisms |2flj | making them hotspots and switching el- 
ements for global chromatin rearrangements. 

We thank R. Bruinsma, K. Kremer, K. Luger, F. 
Miiller-Plathe and J. Widom for helpful discussions. 
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